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ABSTRACT 



We consider a long fiber-optical link consisting of alternating dispersive and nonlinear seg- 
ments, i.e., a split-step model (SSM), in which the dispersion and nonlinearity are completely 
separated. Passage of a soliton (localized pulse) through one cell of the link is described 

Oh . 

by an analytically derived map. Multiple numerical iterations of the map reveal that, at 
values of the system's stepsize (cell's size) L comparable to the pulse's dispersion length zd, 
SSM supports stable propagation of pulses which almost exactly coincide with fundamental 
solitons of the corresponding averaged NLS equation. However, in contrast with the NLS 
equation, the SSM soliton is a strong attractor, i.e., a perturbed soliton rapidly relaxes to 



it, emitting some radiation. A pulse whose initial amplitude is too large splits into two 

O , 

solitons; however, splitting can be suppressed by appropriately chirping the initial pulse. 

•l-H ' 

On the other hand, if the initial amplitude is too small, the pulse turns into a breather, and, 

> 



below a certain threshold, it quickly decays into radiation. If L is essentially larger than zjj , 
the input soliton rapidly rearranges itself into another soliton, with nearly the same area 
but essentially smaller energy. At L still larger, the pulse becomes unstable, with a complex 
system of stability windows found inside the unstable region. Moving solitons are generated 
by lending them a frequency shift, which makes it possible to consider collisions between 
solitons. Except for a case when the phase difference between colliding solitons is < 0.05 ■ n, 
the interaction between them is repulsive. 
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1. Introduction 



A great potential offered by dispersion management schemes for improvement of data trans- 
mission by means of solitons in optical-fiber networks has been well understood (see, e.g., 
a special journal issue! and a collection of articles!). Dispersion management is based on 
periodic compensation of the accumulated dispersion in a nonlinear optical fiber by means 
of additional fiber segments having strong opposite dispersion!. An objective of the present 
work is to introduce and study a related but different scheme (a split-step model, or SSM), 
in which pieces of a dispersive fiber with negligible nonlinearity periodically alternate with 
fiber segments operating close to the zero- dispersion point. In the latter segments, the dis- 
persion is negligible, while nonlinearity is dominating. In the limit when both the dispersive 
and nonlinear segments are very short, the scheme is nothing else but the commonly known 
split-step algorithm for simulations of the nonlinear Schrodinger (NLS) equation, which is 
usually employed in combination with the fast Fourier transform solving the equation at the 
linear stepH; the algorithm can be further modified to adapt it to specific optical models (see, 
e.g., Ref.i). While it is a well established fact that the split-step algorithm yields very accu- 
rate and stable results when simulating solitons, we here aim to consider a model in which 
separation between the dispersive and nonlinear segments is physical, rather than numerical, 
and the lengths of the segments are not small, being comparable to the soliton's dispersion 
length Zd- Understanding possibilities for the existence of split-step solitons (robust solitary 
pulses) and the study of their dynamics (including interactions between them) in a system 
of this type is of a certain fundamental interest. Besides that, the system may also be of 
practical interest for two different reasons. Firstly, the simplest possibility to upgrade a 
linear system so that to make it possible to transmit solitons through it is to periodically 
insert strongly nonlinear fiber segments or similar elements. Secondly, it is necessary to 
understand accuracy limits of the split-step algorithm which can be encountered with the 
increase of the stepsize. A detailed formulation of our model is given below in section 2. 
It is relevant to note that SSM qualitatively resembles a "tandem" model, recently in- 
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troduced by Tornera, which assumes periodic alternation of layers with strong dispersion 
and strong quadratic nonlinearity. As it was demonstrated in Ref.i, this configuration im- 
proves propagation conditions for x solitons in the temporal domain, providing an efficient 
compensation of the group- velocity walkoff. 

We stress that the very existence of stable solitary pulses in a long SSM with a finite 
(nonsmall) stepsize L is not obvious. Although rigorous theoretical considerations of the 
existence problem may be quite difficult, in section 3 we display direct numerical results, 
which show that a very robust pulse can be easily created in SSM. The analysis is based 
on numerical iterations of a map which is produced by the analytical solution of the model 
within one elementary cell. It is found that, in the case L < zd, launching a pulse with 
the initial shape taken as per the fundamental soliton of the usual NLS equation, which is 
obtained by averaging SSM, immediately gives rise to a stationary pulse that appears to 
remain stable for an indefinitely large number of iterations of the split-step map. In the 
cases when the stepsize is essentially larger than the soliton's dispersion length, the initial 
pulse in the form of the NLS soliton emits some radiation and rapidly rearranges itself into 
another soliton with almost the same area but smaller energy. Finally, in the case L 3> Zp, 
the pulses tend to completely decay into radiation. However, in this case we discover a 
system of stability windows on the axis of the parameter L: inside the windows, the soliton 
is still stable. The last window exists at a fairly large value of L. 

In sections 4, we consider the influence of variation of the amplitude of the initial pulse 
and of taking it in a chirped form. If the amplitude is too large, the pulse splits into two 
moving (separating) solitons. However, the splitting can be prevented by chirping the pulse 
(with a right sign of the chirp). The splitting suppression by chirping gives rise to formation 
of a large-amplitude pulse which is quite close to the fundamental soliton. 

In section 5 we study collisions between two identical solitons. Although SSM, unlike 
its NLS counterpart, is not Galilean- invariant, in simulations it is quite easy to generate 
a soliton with a nonzero velocity shift. The interaction strongly depends on the phase 
difference A0 between the solitons. At nearly all the values of A0 they repel each other, 
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while at A0 = the interaction is attractive. A transition from the attraction to repulsion 
occurs at very small values of A0. 

Because the above-mentioned results are obtained numerically, they will be presented in 
the form of various plots and tables illustrating typical situations. It appears that this form 
of the presentation makes it possible to describe the dynamics of pulses in the system quite 
adequately. 

Lastly, we mention that a similar split-step configuration, including only two fiber seg- 
ments, one nonlinear and one dispersive, can be used for various purposes in the form of a 
fiber loop. In particular, it has recently been shown that a loop of this type can operate 
quite efficiently nonlinear mirror for solitonsB. 



2. The model 

The dispersive segment of the fiber is described by the linear Schrodinger equation, 

iu z - ^/3u TT = 0, (1) 

where u, z and r are, respectively, the local amplitude of the electromagnetic waves, propa- 
gation distance, and local time (r = t — z/V, where t and V are the physical time and mean 
group velocity of the wave packet), and (3 is the dispersion coefficient. The main physical 
assumption concerning this part of the system is that its nonlinearity may be neglected 
(for instance, because the fiber's effective cross section is large!). Obviously, Eq. (1) can 
be solved by means of the Fourier transform in r, which is used as a part of the standard 
split-step numerical algorithm!. In an explicit form, the solution is (it can also be obtained 
by means of the Green's function for the linear Schrodinger equation) 



l+isgn/3 /■+«> 
u{z, t) = — . / w(U, rjexp 

!2-k\P\z 



2i t /\2 T / 

-(r-r) dr 



dr'. (2) 



In the nonlinear segment, it is assumed that the dispersion may be neglected (practically, 
the dispersion can be made quite small in dispersion-shifted fiberJl), then one is dealing with 
the dispersionless NLS equation, 



iu z + \u\ 2 u 



(3) 



where the nonlinearity coefficient is normalized to be 1. We also assume that the third-order 
dispersion may be neglected, which is normally true, unless one is going to consider very 
narrow pulsesi. 

The evolution of the field inside the nonlinear segment is described by an obvious solution, 



which is another basic ingredient of the split-step numerical algorithm. We define the sys- 
tem's elementary cell as an interval between midpoints of two neighboring nonlinear seg- 
ments. Thus, the full transformation (map) for the pulse passing one cell can be represented 
as a superposition of two nonlinear phase transformations (§) corresponding to the nonlinear 
half-segments at the cell's edges, and the Fourier transform corresponding to the dispersive 
segment in the middle of the cell. The transformations (0) and (|j) are linked by the condition 
of continuity of u(z, r) at all the junctions between the segments. 

Note that (as well as in the case of the dispersion management) the only dynamical 
invariant of the model is the energy 



A Hamiltonian of the model can be easily defined too, but its value is not conserved at the 
jumps between different segments. 

Models of this type are usually solved with periodic boundary conditions in r. In order 
to prevent radiation (in those cases when it is generated) from re-entering the integration 
domain due to the periodicity, the linear equation (HD was solved directly by means of 
numerical methods in a broad integration interval AT, placing absorbers at its edges (i.e., 
we did not actually use the Green's-function transform (fj)). In fact, the absorbers emulate 
a real physical effect, namely, the fact that the energy emitted with radiation is lost for the 
pulses. 




(4) 




(5) 
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The present model does not take into regard intrinsic fiber losses, amplification, and fil- 
tering. While all these effects are, of course, important in the real communication networks!, 
it is well known that they may be neglected in the first approximation in the case when the 
amplification spacing is smaller than the cell size L (see, e.g., Ref.i, where this issue was 
considered in detail in terms of the dispersion management), although they may become 
important (for instance, giving rise to a specific instability^) when the spacing between the 
amplifiers and filters coincides with L. In any case, the first objective is to study main fea- 
tures of the soliton dynamics in the large-stepsize SSM, while the losses, gain, and filtering 
may be added later. 

If the lengths of the dispersive and nonlinear segments are and L^, the cell size is 
L = Ld + Ltv- Note that the length can be rescaled so that to make \/3\ = 1 in Eq. ([[]). 
Our simulations (many iterations of the above-mentioned map) have shown that in the case 
of the normal dispersion, (3 > 0, no quasi-stable soliton can be found (in complete accord 
with the absence of the bright solitons in the NLS equation with the normal dispersioni) . 
Therefore, in what is following below, we are solely dealing with the case (3 = — 1. 

With regard to the normalizations adopted, an average (distributed) version of the present 
SSM takes an obvious form, 

1 n l2 

iu z + — — — -u TT H — \u \u = 0, (6) 

2{n + 1) n + 1 

where n = L^/Lp is a relative measure of the strength of the nonlinearity in comparison 
with the dispersion. Note, however, that Eqs. (|l|) and @ are separately invariant with 
respect to transformations 

r — > t/Ad,z — ► z/A 2 D , and u — > A^u, z — > z/A 2 N (7) 

with two independent arbitrary scaling factors A^v and Ap. This also transforms the length 
ratio, n — ► (Au/A^n, which may be employed to fix the parameter n. Throughout the 
paper, we fix n = 3. 

In the numerical analysis of SSM, we, will first of all, launch an initial pulse which 
coincides with the fundamental-soliton solution to the average equation @, Uq(t) = 
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(tj I -y/n)sech(^r) , rj being an arbitrary parameter that determines the soliton's amplitude 
and width. Note that, with fixed j3 = — 1 and n — 3, our model, similarly to the usual NLS 
equation, is still invariant against the scale transformation (|7|) with = Ad. Using this 
remaining invariance, we fix rj = 1 in the initial pulse, so that 

uo(t) = (l/\/riJ sechr. (8) 

Besides the initial pulses (||) which correspond to the fundamental-soliton solutions to the 
NLS equation (P), we will also be using initial shapes differing from (§) by taking an arbitrary 
amplitude, and also by adding chirp to the pulse. 

Thus, after fixing all the scales, there remain the single free parameters in the model, 
the stepsize L, to which there may be added parameters of the initial pulse (the chirp and 
relative amplitude) if it is different from the fundamental NLS soliton. 

3. The fundamental soliton in the split-step system 
A. A moderately large stepsize 

We start with the case L = 1, when the stepsize is not essentially different from the dispersion 
length zd ~ (n + 1), defined for the pulse (||) as per the averaged equation (|). Obviously, 
this case is already drastically different from the usual version of SSM with a very small 
stepsize, which is employed in the numerical algorithms for simulations of the NLS equation. 

In Fig. 1 we display a typical numerical solution for this case with the above-mentioned 
fixed value n = 3. As is seen, on a fairly long propagation distance comprising 1500 step- 
sizes (cells) the input (||) generates a very stable pulse propagating virtually without any 
distortion. In the application to optical communications, the physical value of the soliton's 
dispersion length is, typically, ~ 100 kmi, hence, in the present case with n + 1 = 4, the 
interval z = 1500 corresponds to the actual propagation distance ~ 30, 000 km. 

Despite the fact that there is no visible degradation of the pulse in Fig. 1, we stress 
that the existence of the soliton in the model with a finite stepsize in the rigorous sense is 
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not obvious at all, therefore it is necessary to search for a very weak decay of the pulse, 
provided that any decay does take place. To this end, Fig. lb shows the evolution of the 
pulse's energy, E = J slot \u{r)\ 2 dr (cf. Eq. (||)), where the integration is performed over 
a slot which completely covers the pulse but is much smaller than the total size of the 
integration domain, the latter one being AT = 200. As is seen from Fig. lb, after an initial 
transient stage, a systematic loss of energy at an extremely small rate takes place indeed. 
One may interpret this as very weak emission of radiation, plausibly due to the fact that 
the soliton does not exist in the rigorous sense. An extrapolation of the trend seen in Fig. 
lb suggests that the "quasi-soliton" that we are dealing with may persist over an extremely 
long propagation distance corresponding to ~ 5 • 10 6 stepsizes. In the actual simulations 
the soliton kept its shape intact as long as the simulations could be run (up to ~ 5 • 10 4 
stepsizes). 

This and many other runs of the simulations strongly suggest that, in the case L ~ 1, the 
input in the form of a fundamental soliton of the corresponding averaged NLS equation gives 
rise to a very robust stationary pulse that may be regarded as a fundamental soliton of SSM 
(even if the soliton does not exist in the rigorous sense, we may consider it as a virtually 
existing object). The next natural step is to consider an input essentially differing from 
that given by Eq. (|j). There are two straightforward ways to modify the input: changing 
its amplitude by an arbitrary factor A , and/or adding some chirp 6 to it0, so that the 
modified initial pulse is 



It is well known, from both the variational approximation applied to the perturbed soliton 



by radiation losses; eventually, the soliton will assume the static fundamental shape with a 
smaller value of the energyS. However, if the perturbation of the fundamental soliton is too 
strong, it can completely destroy the pulse. In particular, in the case bo = it is known from 




(9) 



in the form (||) and from direct numerical simulations, that in the usual NLS equation the 
input (^) gives rise to long-lived internal vibrations of the soliton, which are slowly damped 



S 



the exact solution obtained by means of the inverse scattering transform that the soliton will 
disappear if A < A t ^ = 1/20 (the variational approximation predicts a higher destruction 



variational approximation predicts that, in the case Aq — 1, the chirp kills the soliton if 
&o — ^thr = V 71 " 2 ' but it is known from numerical simulations that the soliton still survives 
at somewhat larger values of the chirpS. 

Typical results illustrating the evolution of the pulses generated by the unchirped (60 = 0) 
initial shape @ perturbed by A 7^ 1 are displayed in Figs. 2 and 3, in which A takes 
values, respectively, « 1. 73, 0.4^ « 0. 69, and 0.2^ « 0. 346. As is seen in Fig. 
2a, in the case Aq > 1 the pulse quickly rearranges itself, after several vibrations, into a 
quasistationary state, which can be verified to be very close to the fundamental soliton. The 
rearrangement is accompanied by a considerable loss of energy, which is radiated away. In 
fact, the radiation losses make the fundamental soliton in SSM an effective attractor, similar 
to those in dissipative systems, although there are no direct losses in the present model. 

The strength of the soliton in SSM as the effective attractor is clearly seen in comparison 
with its NLS counterpart. To this end, in Fig. 2c and 2d we display results of direct 
simulations of the averaged NLS equation (|5|) with exactly the same initial pulse as in 
Fig. 2a (in fact, the integration of the NLS equation was also performed by means of 
the split-step algorithm, but with a very small stepsize). A drastic difference between the 
relaxation of the deformed soliton in two models is obvious, despite the fact that the shapes 
of the fundamental soliton in these models are virtually identical. A natural explanation 
to this is the fact that SSM with L ~ 1 is strongly nonintegrable, on the contrary to the 
integrable NLS equation, and, as it is well known, irreversible dynamical processes (such 
as the separation of a soliton and radiation) are much faster in nonintegrable systems^. 
The fast relaxation of the perturbed soliton in SSM has an apparent advantage for the 
applications to telecommunications; on the other hand, fast relaxation assumes a burst of 
radiation emission. In a more realistic model, including optical filters, the radiation will be 
absorbed by the filters. 




When the chirp is too strong, it also destroys the solitontll. The 
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A qualitatively different result (splitting of the initial pulse into two) is produced by 
the multiplication of the input by a still larger factor A . This effect will be considered 
separately in the next section. 

The perturbation of the initial pulse by the multiplier Aq which is considerably smaller 
than 1 produces a very different effect, transforming the stationary soliton into a breather 
with long-period undamped internal vibrations (Fig. 3). Finally, if Aq is too small, this 
results in quick destruction of the pulse (Fig. 4). Thus, a threshold value of A separating 
the cases when the soliton survives and those when it perishes is located somewhere between 
0.69 and 0.346 (recall that A^ = 1/2 for the NLS equation). An accurate value of A t hr can 
be found in a straightforward way by means of longer simulations. 

A different generic type of the perturbation is chirping the input. A typical example 
with a strong chirp is shown in Fig. 4. In fact, this is another manifestation of the remark- 
able robustness of the soliton in SSM: the pulse quickly sheds off considerable amounts of 
radiation, loosing with it nearly half of its initial energy (Fig. 4b), and again rearranges 
itself into a quasistationary state. In Fig. 4c, we compare the shapes of the central parts 
of the input and output pulses, along with the chirp distributions in them (the local chirp 
shown in Fig. 4c is TT , where 4>(z,t) is the internal phase of the pulse). Note that the 
final pulse has the amplitude and width which are simultaneously smaller than those of the 
initial pulse (this is possible in view of the considerable radiative losses). As for the chirp in 
the final pulse, two features are noteworthy: its absolute values are much smaller than the 
chirp 2b = 2 in the initial pulse, and, moreover, the average value of the chirp across the 
central part of the output pulse is nearly zero. Thus, we conclude that SSM quite efficiently 
suppresses the chirp, which is not the case in the usual NLS equation. In fact, because the 
chirp is a dynamical variable conjugate to the soliton's width, within the framework of the 
variational approximationS, the chirp suppression is closely related to the above-mentioned 
fast suppression of the pulse's internal vibrations in SSM (Fig. 2a). Nevertheless, simula- 
tions with very large values of bo (not shown here) demonstrate that a very strong chirp 
destroys the pulse, as one would expect from the analogy with the NLS equation. 
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B. A very large stepsize 



In the case L>1, i.e., when the stepsize is considerably larger than the soliton's dispersion 
length, the situation is very different from that described above. This case will take place 
if, in order to increase the data transmission rate, narrower solitons are launched into the 
fiber linkSi (note that z D scales inversely proportional to the soliton's temporal widthi). 

A drastic difference revealed by the numerical computations in the case L 3> 1 is that 
the initial pulse in the form of the NLS fundamental soliton (|8|) no longer remains a quasis- 
tationary solution of SSM. Instead, as it is illustrated by Fig. 5 pertaining to L = 10, the 
initial pulse undergoes fast evolution, decreasing its amplitude and getting broader. Never- 
theless, this does not lead to decay of the pulse, and it keeps more than half of its initial 
energy (Fig. 5b) . The outcome pulse can be very accurately fitted to the soliton-like ansatz, 



with some amplitude A and width T, and Fig. 5d shows that the pulse's chirp (taken at a 
value of z corresponding to the junction between two cells) remains very small. 

The same ansatz (|10[) provides for a good fit for the outcome solitons generated by 
simulations at other large values of L (at which the solitons are still stable, see below). The 
values of the fitting parameters for several large stepsizes L, all corresponding to n = 3 and 
the initial pulse in the form @, are collected in Table 1. In the same table, we also give 
values of the numerically measured area of the pulse corresponding to the expression (|10|), 



(unlike the energy (|5|), the area is not a model's dynamical invariant). A well-known theorem 
states that an arbitrary pulse which is considered as an initial condition to the averaged NLS 



u\ = A sech (t/T) , 



(10) 




(11) 



equation (|6|) cannot give rise to a soliton in this equation unless the pulse's area exceeds 
a threshold valueEl 5 m i n = rr x l 2 In ^2 + \^3j ~ 1-317 • n -1 / 2 ; in particular, the area of the 
NLS fundamental soliton (|8|) is S'soi = 7m -1 / 2 . 
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It is clearly seen from Table 1 that the soliton's area remains fairly close to (slightly larger 
than) the area of the NLS soliton, S so \ = n/y/3 ~ I. 81 (recall we have set n = 3). Thus, we 
infer that launching a fundamental NLS soliton into SSM with a large stepsize L transforms 
the soliton (unless it is unstable, see below) into another one which belongs to the same 
family of the fundamental solitons, having nearly the same area, but with an essentially 
smaller energy. The data from Table 1 suggests that, roughly, the residual energy of the 
outcome fundamental soliton decreases ~ 1/ y/L, although there are considerable deviations 
from this power law. Such a transformation of one fundamental soliton into the other is a 
manifestation of strong nonintegrability of SSM with large L. 

It may happen that the outcome pulse obtained at large L, i.e., a fundamental soliton 
of SSM proper, does have a systematic difference in its shape from its NLS counterpart. 
For instance, we tried to fit the data to an ansatz more sophisticated than (0), viz., \u\ = 
A [sech (r/T)] Q . The result was that a may take values between 1 and 1.6, having then the 
width somewhat larger than that given in Table 1 . However, the available accuracy of the 
numerical data does not allow us to make a decisive conclusion that a > 1 provides for a 
really better fit than the simplest ansatz ([TDD . 

Lastly, we stress that it may be quite plausible that SSM does not have any soliton 
solution in the rigorous sense. However, in practical terms, it is quite easy to distinguish 
between cases when virtually stable solitons persist on extremely long propagation distances, 
and those when any pulse quickly decays, as it is shown in detail in the next subsection. 

C. Stability windows 

As it was already shown above, a general trend is that, with the increase of L, the amplitude 
of the persisting soliton decreases, and at extremely large values of L solitons decay into 
radiation. However, detailed simulations reveal a very nontrivial feature: there are stability 
windows on the axis of the parameter L, alternating with regions in which the soliton is 
either completely unstable or "semi-unstable" . The latter means that the soliton suddenly 
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jumps down to a smaller amplitude, at which it seems to persist. Increasing L by steps 
AL = 1, we have found that the solitons remain continuously stable at L < 14; then, the 
instability or semi-instability takes place at 



Accordingly, stability windows were found at 

L = 18, L = 20, L = 22, L = 24, L = 26, L = 36, L = 51, andL = 59 . 

To illustrate these features, we display several examples in Fig. 6: the last representative 
of the continuously stable solitons at L = 14 (a), the first semi-unstable one at L = 15 (b), 
the first fully unstable case at L = 16 (c), and the first stability window at L = 18 (d). We 
do not show here solitons belonging to the farthest stability windows, as it is rather difficult 
to display them using the same scale as in Fig. 6. However, in Fig. 7 we separately show 
the eventual shape of the soliton in the last stability window (L = 59). 

The above system of the stability windows was found with a rather crude resolution, 
AL = 1; we expect that using a finer resolution would reveal a much more complex system 
of windows, and it seems plausible that their total number is infinitely large, i.e., the window 
system may have a fractal structure. 

Lastly, we notice that, in the stability windows corresponding to the largest values of L, 
the soliton does not get completely separated from radiation, which is illustrated by Fig. 7 
(cf. Fig. 5c pertaining to L = 10, where there is no visible overlapping with radiation). 

A complex system of alternating windows was earlier discovered in a very different non- 
linear model based on the Goldstone (also called 4 ) equation for a real function (p(x,t), 



Numerical simulations of collisions between topological solitons (kinks) with opposite po- 
larities and opposite velocities ±v in Eq. ( |13|) had revealed that the collision results in 



15 < L < 17, L = 19, L = 21, L = 23, L = 25, 27 < L < 34, 



37 < L < 50, 52 < L < 58, and L > 60 . 



(12) 



4>tt - 4>xx - 4> + 3 = o . 



(13) 
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annihilation of the kinks into a breather (which is subject to subsequent slow radiative de- 
cay) if the velocity is very small, and the collision is quasi-elastic if v is sufficiently close to 
the limiting velocity t> max = 1. Between these two system of alternating windows 

for annihilation and quasi-elastic collisions was found (see Ref.i and references therein). 
A semi-quantitative explanation to these findings was based on consideration of exchange 
between the kink's kinetic energy and the energy absorbed by an internal oscillatory degree 
of freedom which the kink is known to haveS. It may happen that the SSM soliton also has 
an internal degree of freedom, which may be related to the system of the stability/instability 
windows. However, a detailed explanation of this feature should be a subject of a separate 
work. 

To conclude this section, it may be also relevant to mention that, in the studies of the 
dispersion-management models, it has been recently found that stable propagation of a 
soliton is possible when the average dispersion takes a slightly normal value (corresponding 
to (3 > in Eq. (|l]))i'i@. On the basis of additional simulations, we have concluded that, 
in the present model, no (quasi) stable soliton can be found if the dispersion is normal. 



4. Splitting of a large initial pulse and its suppression by means of chirping 

Getting back to the "moderate case" with L — 1, we now aim to consider the evolution 
of the pulse (|9|) with large values of the amplitude A Q . In this case, the most interesting 
observation is spontaneous splitting of the pulse into two "splinters" when Aq exceeds a 
certain threshold value, which is close to 2 if L — 1. A typical example of the onset of 
the splitting regime is shown in Fig. 8. In the figure, we keep A Q = 2.2, but increase L 
from 1/3 to 1. Note that in the NLS equation, which corresponds to the limit L — > 0, the 
unchirped {b = 0) pulse @ does not split at arbitrarily large values of A ; instead, it gives 
rise to a higher-order soliton (breather )0. In contrast with this, Fig. 8 clearly suggests that, 
with the increase of L, i.e., with increasing departure from the NLS limit, the SSM pulse 
with a sufficiently large amplitude and without initial chirp performs a transition from a 
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strongly pulsating breather to splitting. The splitting in the absence of chirp seems to be a 
characteristic feature of nonintegrable models, as it is definitely impossible in the integrable 
NLS equation. Note that a similar feature was recently observed in simulations of a different 
nonintegrable model, viz., the cubic-quintic NLS equation^. 

It is well known that initial chirp can easily split the pulse even in the integrable NLS 
equation (see, e.g., Ref.0). Therefore, it is natural to assume that adding the chirp (with a 
proper sign) to the initial pulse which splits in the absence of chirp may balance the trend 
to the splitting and, eventually, suppress it, providing for the generation of a quasistationary 
nearly fundamental soliton. Simulations show that this may indeed be the case, a typical 
example of which is displayed in Fig. 9. As is seen in Fig. 9b, more than half of the initial 
energy of the pulse is lost with emitted radiation, but one may say that a single soliton is 
eventually formed, with a small "satellite" attached to it, cf. Fig. 8c. Note that the average 
value of the chirp in the central body of the output pulse (Fig. 9c) is close to zero. Thus, 
properly selected prechirping of the large-amplitude pulse not only can prevent its splitting, 
but also helps to transform it into a nearly fundamental soliton, rather than leaving it in a 
strongly vibrating state (cf. Figs. 9a and 8a,b). 

On the other hand, adding the chirp with the opposite sign to the initial large-amplitude 
pulse, we observed that the splitting, instead of being suppressed, is strongly enhanced (not 
shown here). It is, of course, quite natural that the effect of chirping crucially depends on 
its sign. 

5. Collisions between solitons in the split-step system 

Equation (0) and, hence, SSM as a whole, is not Galilean invariant, unlike the classical 
NLS equation. Nevertheless, a soliton with a nonzero shift of its (inverse) velocity can be 
generated by means of a simple numerical trick0: a known solution for the "quiescent" 
soliton should be multiplied, at the initial point z = 0, by exp(— iujt) with an arbitrary 
frequency shift u, which can generate a velocity shift via the dispersive part of the model. 
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Thus, we can easily generate two solitons with opposite velocities and simulate their collision 
(in fact, moving solitons generated by the splitting of a large-amplitude initial pulse were 
already observed above in Fig. 8). 

If the frequency shift u is large enough, the collision seems quite elastic, see a typical 
example in Fig. 10. At smaller frequency shifts, however, collisions may give rise to very 
different results. Below, we display a set of typical examples obtained with u = ±0.05, 
while L — 1 and n = 3 (the same values for which a majority of the above results have been 
presented) . 

The outcome of the collision in the case of small u strongly depends on the phase differ- 
ence A0 between the colliding solitons (we consider only collisions between identical ones). 
As is well known, the interaction between solitons of the NLS type is attractive if their phase 
difference is zero, and repulsive if A0 = ir. In accord with this, we observed that the solitons 
readily pass through each other - obviously, due to the attraction - in the case A0 = (Fig. 
11a), while in the cases A0 = tt and A0 = ir/2 they clearly repel each other and avoid 
mutual passage (see Figs, lib and 11c). More detailed simulations have demonstrated that, 
in fact, the interaction between the solitons is attractive only if the phase difference between 
them is very small, < 0.05 • ir. 

6. Conclusion 

We have introduced a model of a long fiber-optical link consisting of alternating dispersive 
and nonlinear segments, i.e., a model with a complete separation of the dispersion and 
nonlinearity. Losses and gain were assumed to be compensated, therefore they were not 
included into the model. 

Passage of a soliton (localized pulse) through one cell of the system was described by 
an analytically derived map. Multiple numerical iterations of the map reveal that, at val- 
ues of the system's stepsize L comparable to the soliton's dispersion length zd, the system 
supports indefinitely long stable propagation of pulses which almost exactly coincide with 
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fundamental solitons of the corresponding averaged NLS equation. If the initial pulse is 
perturbed, it quickly relaxes into the fundamental state, shedding off some radiation. The 
relaxation is much faster than in the averaged NLS equation, which allows one to consider 
the soliton in the present model as an effective attractor. This property has an advantage 
for telecommunications. However, an initial pulse whose amplitude is too large splits into 
two moving solitons. The splitting could be suppressed by appropriately chirping the ini- 
tial pulse, so that a nearly fundamental soliton appears again. On the other hand, if the 
amplitude of the input pulse is too small, it turns into a long-period breather, and, below a 
certain threshold, it quickly decays into radiation. 

If L is essentially larger than zd, an input pulse in the form of an NLS soliton rapidly 
rearranges itself into another soliton, whose area remains nearly the same as that of the 
NLS fundamental soliton, but the energy is considerably smaller than that of the input 
pulse. With further increase of L, the pulse becomes unstable; however, a complex system 
of stability windows was found inside the unstable region. 

Moving solitons could be prepared by initially giving them a frequency shift, which makes 
it possible to simulate collisions between solitons. A result is that, except for a case when 
the phase difference A0 between colliding identical solitons is very small, the interaction 
between the solitons is repulsive, and they avoid passing through each other, unless their 
relative velocity is very large. 
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L 


5 10 14 20 26 36 51 59 


A 


0.45 0.33 0.28 0.24 0.20 0.16 0.12 0.10 


T 


1.28 1.88 2.15 2.50 2.98 3.60 5.10 5.90 


S 


1.83 1.87 1.87 1.85 1.87 1.83 1.84 1.84 



Table 1. Values of the fitting parameters from the ansatz ( ]T0| ) for the outcome soliton 
in the split-step system with n = 3 and large values of the stepsize L (the initial pulse was 
always taken in the form (|8|)). Also given in the table are values of the soliton's area ([TTJ) 
(note that, for n — 3, the area of the NLS soliton is S so \ = n/VS ~ 1. 81). 
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Figure captions 

Fig. 1. The numerical solution of the split-step model with the stepsize L = 1 and the 
ratio of the lengths of the nonlinear and dispersion segments n = 3. The initial pulse is the 
fundamental soliton @ of the corresponding averaged NLS equation (|]). Shown are (a) \u\ 2 
vs. z and r, and (b) the evolution of the pulse's energy over 1500 stepsizes. 

Fig. 2. The same as in Fig. 1 in the case when the perturbed initial pulse ( [5J) is taken 
with bo = and Aq = 1.73; the panels (c) and (d) additionally show evolution of the same 
pulse and its energy in the averaged NLS counterpart (|J) of the split-step model. 

Fig. 3. The same as in Fig. la in the case when the initial pulse is multiplied by 0.69 
(a) or 0.346 (b). 

Fig. 4. The same as in Fig. 1 in the case when the chirp, bo — 1, is added to the initial 
pulse (H) with Ao — 1. The panel (c) additionally shows the initial (dotted) and final (solid) 
shapes of the central part of the pulse, |w(t)| 2 , and the distribution of chirp, therein. 

Fig. 5. The same as in Fig. 1 with L = 10 and n = 3. The panel (c) additionally 
shows the best fit of the shape of the final pulse \u(r)\ to the ansatz (|T0|), and (d) is the 
distribution of chirp in the central part of the output pulse. 

Fig. 6. Alternation of stable and unstable pulses at large values of the stepsize (for 
n = 3): (a) L = 14 (the end of the continuous stability region); (b) L = 15 (the first 
semi-unstable case); (c) L = 16 (the first fully unstable case); (d) L = 18 (the first isolated 
stability window). 

Fig. 7. The eventual shape of the fundamental soliton, |«(r)|, at L = 59 (with n — 3). 
The dotted curve is the best fit to the ansatz (|i~0|). 

Fig. 8. Splitting of the pulse generated by the initial shape @ with bo = and large 
additional multiplier A = 2.2 in the model with n = 3 and different values of the stepsize: 
(a) L = 1/3; (b) L = 1/2; (c) L = 1. 

Fig. 9. Suppression of the pulse splitting shown in Fig. 8 by prechirping (with b = 1) 
the initial shape (|9|): (a) \u\ 2 vs. z and r; (b) the evolution of the pulse's energy; (c) the 
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distribution of the chirp in the output field (taken at z — 300), 

Fig. 10. A quasi-elastic collision of two identical solitons in the split-step model with 
L — 1, n — 3, and the phase difference A<p = 0. The solitons were set in motion by giving 
them large initial frequency shifts uj = ±0.76. 

Fig. 11. Collisions between two identical solitons that were "pushed" by small initial 
frequency shifts u = ±0.05 in the case L — 1 and n — 3. The phase difference between the 
solitons is A0 = (a), A0 = 7r (b), and A0 = 7r/2 (c). Because in the case of the small 
relative velocity it is difficult to show collisions by means of three-dimensional plots similar 
to that shown in Fig. 10, we here instead display contour plots, that illustrate the collisions 
process quite clearly. 
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